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Striped superconductors in the extended Hubbard model. 



O 

o 
o 

(N 
< 

OO 



Ivar Martin, Gerardo Ortiz, A. V. Balatsky, and A. R. Bishop 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 

(February 1, 2008) 

We present a minimal model of a doped Mott insulator that simultaneously supports antiferro- 
magnetic stripes and d-wave superconductivity. We explore the implications for the global phase 
diagram of the superconducting cuprates. At the unrestricted mean-field level, the various phases of 
the cuprates, including weak and strong pseudogap phases, and two different types of superconduc- 
tivity in the underdoped and the overdoped regimes, find a natural interpretation. We argue that on 
the underdoped side, the superconductor is intrinsically inhomogeneous - striped coexistence of of 
superconductivity and magnetism - and global phase coherence is achieved through Josephson-like 
coupling of the superconducting stripes. On the overdoped side, the state is overall homogeneous 
and the superconductivity is of the classical BCS type. 
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Experimental evidence increasingly suggests that mi- 
croscopic inhomogeneous "stripe" states are ubiquitous 
in the doped cuprates, as well as in other complex elec- 
tronic materials B. Nanoscale stripe morphologies have 
been inferred in YBCO and LSCO from neutron scat- 
tering and angle resolved photoemission experiments. In 
the superconducting phase, the stripes appear to coexist 
with superconductivity in a range of dopings without de- 
stroying the global phase coherence. The main issue that 
we address in this Letter is the nature of the supercon- 
ducting state in the presence of stripes. 

Inhomogeneous quantum states in non-superconducting 
lattice models are no less common than they are in the 
experiments. Many lattice models possessing antiferro- 
magnetic (AF) ground states at half-filling, when doped 
away from that filling develop stripes at the unrestricted 
mean- field (MF) level ||-||. Exact solutions may lead 
to fluctuations that introduce dynamics into the MF 
solutions, but are expected to preserve the qualitative 
features. The lower the spatial dimension the more im- 
portant the quantum fluctuations are and, sometimes, 
MF solutions do not reproduce the exact large distance 
physics of the problem. However, often they pick up the 
low lying manifold of excited states which becomes rel- 
evant at low enough temperatures. Also, 3-dimensional 
coupling in real materials helps to reduce the effect of 
fluctuations. 

We consider here a minimal model with stripes to il- 
lustrate our conclusions. We employ the 2-dimensional 
one-band Hubbard Hamiltonian with hopping t and on- 
site repulsion U Superconductivity is introduced by 
including the nearest neighbor attraction V , which pro- 
duces pairing predominantly in the d-wave channel close 
to half- filling ||. The effective minimal Hamiltonian is 
thus 

H = -t c\* c j<7 + u ^n^nn + V^n^-, (1) 

' (ij) 

where the operator c\ a [cj a ) creates (annihilates) an 
electron with spin a on the lattice site i, and n\ — 



C If C if 



c ll c il represents the electron density on site i. 



For our computations, we use the unrestricted mean-field 
approximation to this Hamiltonian, 



H MF = ~* X] c i" c j° + U Y1 n ^( n n) + ( n it) n H 

(ij)a i 



H.c. 



(2) 



where Ay = V(c,-;Cjf) is the MF superconducting order 
parameter. The direct Hartree terms in V are neglected 
since the magnitude of the effective nearest neighbor at- 
traction is expected to be much smaller than the on-site 
repulsion U. Hence, it should not affect the diagonal part 
of the Hamiltonian, which is responsible for the charge 
and spin order. Therefore, the effect of V in our model 
is limited to the generation of superconducting correla- 
tions. We do not address the very important issue of the 
microscopic origin of the attraction V. Our goal is only 
to construct a minimal model that may help to qualita- 
tively understand the rich phase diagram of the cuprates. 
Unless stated otherwise, standard parameter values used 
are U — At and V = — Q.9t. This choice of parameters al- 
lows us to clearly demonstrate our conclusions regarding 
the inhomogeneous superconducting phase. 

To self-consistently solve the MF equations we use 
an iterative scheme. This consists of two stages which 
are repeated until convergence is achieved: (1) Diago- 
nalization of the MF Hamiltonian, and (2) Update of 
the MF parameters. An important feature of our ap- 
proach is that all physical quantities are allowed to vary 
from one lattice site to another, e.g., (rijf) ^ (nil) and 
{n%o) ^ (ni+aa)- Generically, we are seeking inhomo- 
geneous solutions whose typical correlation lengths £ in- 
volve several lattice spacings. Therefore, it is important 
that the simulated supercell size N x x N y (with peri- 
odic boundary conditions) is such that N x ^ y > £ x>v . The 
Hamiltonian in Eq. (||) can be rewritten in the matrix 
form, Hmf = He, with c = (c, t , e|^, c^, cj^, • • • ) T , 
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where H is a (2N x N y ) x (2N x N y ) hermitian matrix. By 
applying a unitary transformation a (a -1 = a*) the 
Hamiltonian matrix can be diagonalized as H — aDa -1 , 
with D nm = 5 nm E n . The Hamiltonian can be diagonal- 
ized in the Bogoliubov quasiparticles, 



7n 



a nm C n 



(3) 



with energies E n . By reexpressing the original creation- 
annihilation operators in terms of the Bogoliubov quasi- 
particles, one can recompute the parameters of the MF 
Hamiltonian. For example, 

(n iT ) = (cj T C iT ) = J2( a h,nl f n a th™lrn} 
nm 

= \a ihn \ 2 n F (E n ), (4) 

n 

where np(E n ) is the Fermi-Dirac distribution function. 
Repeated until the convergence, the iterations produce 
the spatial profiles of the self-consistent density and or- 
der parameter. 

A typical zero-temperature MF inhomogeneous solu- 
tion is shown in Fig. [l]. In the lowest energy configu- 
ration, the spin density develops a soliton-like AF anti- 
phase domain boundary — a stripe — at which the AF 
order parameter changes sign. At the domain boundary, 
the electronic charge density is depleted. The width of 
the domain wall, £,DWi decreases with increasing on-site 
repulsion U . However, for values of U that are not much 
larger than the hopping t, the charge per unit length of 
the optimal (the lowest energy) stripe remains the same 
and is close to unity near half- filling. This result, first 
demonstrated in this simple model by Schultz 0, is a 
direct consequence of doping-dependent nesting in the 
Hubbard model. The bond-centered stripes are favored 
relative to the site-centered ones, although the energy 
difference in our case is small due to the smooth charge 
distribution. For a different band structure the exact 
relation between the doping x and inter-stripe distance, 
L{x), may change; however, any model whose ground 
state is AF at zero doping, is expected to have AF stripes 
for a finite doping, with incommensuration proportional 
to the doping, X/L(x) cx x, near half-filling. For ex- 
ample, negative next-nearest neighbor hopping t' (rele- 
vant in the hole-doped cuprates M), modifies the stripe 
filling without compromising the stripe phase stability 
relative to commensurate AF at the MF level |Q. The 
stripe filling is a monotonically decreasing function of 
the magnitude of t', with the magic filling 1/2 occurring 
when t' — — 0.35t. Experimentally, however, the value of 
the effective doping-dependent t' is still not well defined. 
Among the consciences of the doping dependence of t' 
are doping-dependent stripe filling, and a possibility of a 
transition from vertical to diagonal stripes as a function 
of doping. While we focus here on the model case t' = 0, 



our conclusions remain qualitatively the same also in the 
presence of a moderate negative \t'\ < 0.5t. 

For large enough doping levels, such that L(x) < £dw, 
the AF stripes begin to overlap. In this regime the ex- 
citation spectrum is no longer fully gapped and mobile 
carriers appear. Further doping mostly changes the am- 
plitude of the spin and charge density waves, only slowly 
modifying the incommensuration, 1/L(x) [§). When the 
stripes are sufficiently close to melting, the AF aspect 
of the problem becomes unimportant, and the supercon- 
ductivity is of the conventional BCS type. 
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FIG. 1. Typical example of density and superconducting 
order parameter profiles in a stripe state (here period 17). 
The top two bar charts represent the site-dependent spin 
and charge densities, respectively. The contour plots indi- 
cate the sites with low (blue) and high (red) values of the 
corresponding densities. The bottom four plots show the 
values of the superconducting order parameters, defined as 
Af (s ' = (A iir i 9ht +A iiie / t =FAi, tl p^Ai idotun )/4 for d-wave (ex- 
tended s-wave) order parameter on site i (U = At, V = — 0.9i). 
Different choices of parameters lead to qualitatively similar 
patterns, with stronger U leading to a stronger AF order and 
more attractive V causing the superconducting stripes to be- 
come wider and larger in amplitude. The doping level is 5.9%. 



The superconducting order parameter A^- s ^ is max- 
imal on the stripes and is not smooth (even within the 
stripe) due to the presence of the AF background. This 
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happens since the order parameter is sensitive to the spin 
density on sites i and j. If i belongs to the spin-down 
sublattice and the neighbor j is on the spin-up sublat- 
tice, then the order parameter is large, and vice versa. 
Notice, that in addition to the dominant <i-wave compo- 
nent, there is a small extended s-wave component gener- 
ated on the stripe, which can be interpreted as a distor- 
tion of the d-wave at the level of about 10%. This hap- 
pens because a symmetry of the lattice has been broken. 
The superconducting stripe is pinned by the AF phase 
domain boundary. Since the spin density on the domain 
boundary is small, a natural interpretation is that the 
superconductivity is suppressed in the regions of large 
AF order. Nevertheless, the width of a superconduct- 
ing stripe, £sc, is determined not only by the width of 
the AF stripe, but also increases rapidly with increasing 
magnitude of the nearest neighbor attraction, V. For our 
choice of parameters, ^dw ~ 4 and ^sc ~ 8 lattice sites. 
For dopings smaller than about 10% (corresponding to 
L(x) > 10 lattice sites) the stripes have negligible over- 
lap. In this regime, the amplitude of the superconducting 
order parameter on the stripes no longer depends upon 
the stripe-stripe separation. For higher doping levels, an 
overlap between the superconducting order parameters 
on adjacent stripes is established. 



A central question is connected to the conducting prop- 
erties of the resulting inhomogeneous state: Is it a global 
superconductor, a metal, an insulator, or some unusual 
anisotropic phase? To resolve this issue we use the con- 
cepts of charge stiffness D c Q , which measures the sensi- 
tivity of the ground state to changes in boundary condi- 
tions, and anomalous flux quantization [ jTo[ , which pro- 
vides a direct signature of the Meissner effect. These 
concepts, together with topological quantum numbers 
|pr| , are routinely employed to study the localization 
properties of models of interacting electrons. To com- 
pute D c one needs to determine how the energy of a 
system with a fixed number of particles, E, depends on 
the twist in the boundary conditions, 9 E [0, 2ir). The 
twist of the boundary conditions is independently ap- 
plied along each spatial direction j = x,y and implies 
that CNj+i = exp(i8)ci. The special case of = cor- 
responds to strictly periodic boundary conditions. Text- 
book schematics of the energy dependence, E(Q), are 
shown in Fig. ||A. The calculated many-particle spec- 
trum E(Q) for a system with a stripe separation of 17 
lattice sites in our model is shown in Fig. |^B. The energy 
curves imply that the system is superconducting, however 
the superfluid stiffnesses along and across the stripes are 
drastically different. However, for a smaller stripe peri- 
odicity (Fig. [2jC) , due to substantial overlap between the 
stripes, superconductivity is almost as strong across the 
stripes as it is along the stripes. 
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FIG. 2. (A) Schematics of typical energy spectra, E(Q). 
Insulating behavior is characterized by a smooth curve of pe- 
riod 27r, with an amplitude that decays exponentially with 
system size. The metallic and the superconducting states 
typically have patterns with cusps, with variations in E(Q) 
that decrease as a power of the system size in such a way 
that the charge stiffness, D c — j^d 2 E/dQ 2 , remains con- 
stant. The superconducting behavior is distinguished from 
the metallic one by the reduced periodicity, referred to as 
anomalous flux quantization (AFQ). There is a direct corre- 
spondence between AFQ and the fact that the minimum flux 
that can penetrate a superconducting system is one half of 
the flux quantum, <f>o/2 = hc/2e. The exact period it is only 
achieved in the thermodynamic limit; for a finite system, there 
is only a signature of the reduced period. (B) Calculated en- 
ergy spectrum E(Q) for a system of size 10 x 17 (see Fig. |l|). 
Along the stripes there is a pronounced AFQ signature, which 
implies that the system is superconducting along the stripes. 
The energies were computed at a fixed chemical potential. To 
convert to the energy at a fixed average number of particles, a 
density adjustment by the amount —fj,(n(Q) — n(0)) has been 
made. Across the stripes, the stiffness is very small, but has 
a period of 7r. Hence the system appears to be a global super- 
conductor, but with an extremely small superfluid stiffness in 
the direction across the stripes. (C) Calculated energy spec- 
trum E(Q) for a system of size 10 x 15 with 3 stripes of period 
5. Due to substantial overlap between the superconducting 
stripes, this system has superconducting strengths that are 
comparable in both directions. 



Arrays of superconductors separated by insulating re- 
gions, known as Josephson junction arrays, have non- 
trivial conducting properties. Depending on the rela- 
tive strength of the coupling between the superconduc- 
tors and the charging energy, such systems can be ei- 
ther superconductors or insulators At the MF level, 
the charging aspect of the problem is absent, and hence 
we can expect a global superconductor at zero tempera- 
ture no matter how weak the interchain coupling is. In 
a real system, for a sufficiently weak coupling the su- 
perconducting behavior will be suppressed. In the case 
of superconducting channels separated by AF insulators, 
the coupling is inversely related to L{x). For very large 



3 



LA-UR-OO-XXXX 



February 1, 2008 



submitted to Nature 



L(x) at low doping, L(x) 3> £sc, the overlap between the 
superconducting stripes, and hence the superconducting 
transition temperature T c , is exponentially small. As the 
distance between the stripes decreases (larger doping), 
the overlap of the superconducting condensate wave func- 
tions should establish a phase coherent superconducting 
state. Indeed, this is qualitatively what we observe al- 
ready at the mean-field level in the striped supercon- 
ductors. For a large superconducting stripe overlap, the 
effective coupling between the stripes is non-exponential. 
In this regime, the experimentally measured supercon- 
ducting transition temperature is proportional to incom- 
mensuration |]l3| , which implies that the effective Joseph- 
son coupling scales as 1/L(x). 

A possible experimental test of the Josephson-coupled 
superconductor scenario proposed here (see also flij ) can 
be performed by measuring the in-plane Joscphson plas- 
mon resonance. The resonance should be present in the 
microwave-frequency range and is excitable by an in- 
planc electric field. 
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FIG. 3. Schematic phase diagram obtained by super- 
imposing the antiferromagnetic (AF) / striped (ICAF) 
and the d-wave superconducting (SC) phase diagrams. 
In the intersection region we distinguish the subregions 
of Josephson-coupled striped superconductor (SSC), and 
non-superconducting "strange metal" (SM), which is neither 
a superconductor, nor a simple insulator. The upper bound- 
ary of the AF/ICAF corresponds to the weak pseudogap 
crossover, and the line between the pure AF/ICAF and the 
SM marks the strong pseudogap crossover. In the very low 
doping region, the superconducting aspect of the problem be- 
comes irrelevant, and the phase diagram is dominated by the 
physics of antiferromagnets. A detailed finite-temperature 
study is required to precisely locate the left boundary of the 
SSC region, and hence to determine if the critical point P is 
indeed a penta-critical point. 



From our zero-temperature analysis of the coexistence 
of AF stripes (ICAF) and superconductivity, a simple 
qualitative thermodynamic phase diagram emerges. In 
the conjectured phase diagram, we utilize the finite- 
temperature AF/ICAF phase diagram of the Hubbard 
model together with the superconducting (SC) phase di- 



agram of the t-V model. The SC phase diagram is ob- 
tained in the homogeneous MF 0, while the AF/ICAF 
phase boundary is constructed under the assumption of 
the second order phase transition between the homoge- 
neous and inhomogeneous states || . For a suitable choice 
of parameters, for instance U = 2t and V = — t, the SC 
and the AF/ICAF regions in the phase diagram inter- 
sect, as shown in Fig. || The boundary between the 
AF and ICAF phases corresponds to an infinite period 
stripe modulation, and implies that the incommensura- 
tion is a decreasing function of temperature. The energy 
scale associated with the AF/ICAF region of the phase 
diagram is much larger than that of the SC part. Thus, 
one expects that only the SC phase boundary is modified 
when it passes through the AF/ICAF region. The cen- 
tral result of our work is that the superconductivity does 
not disappear in the region of the AF stripes, but rather 
becomes striped, with anisotropic superfluid stiffness. 

Based on familiar Josephson coupling physics, in the 
region of coexistence of superconductivity and stripes, we 
can expect a part that is a globally coherent striped su- 
perconductor (SSC). The rest of the intersection region is 
covered by an exotic phase which, if it were perfectly ori- 
entationally ordered, would be a superconductor in one 
direction and a strongly-correlated insulator in the other. 
In reality, due to the meandering of the stripes and their 
break-up into finite segments |H| , the state is likely to be 
highly inhomogeneous and neither an insulator, nor a su- 
perconductor, but also not a simple metal. In agreement 
with the experimental attribution, we refer to this region 
as a "strange metal" (SM). The line separating the SM 
from the AF/ICAF region, in the context of the experi- 
ments, can be associated with the crossover to the strong 
pseudogap regime, and corresponds to the opening of the 
superconducting gap. The high-temperature AF/ICAF 
phase boundary marks to onset of the weak pseudogap. 
For very small dopings the MF stripe separation becomes 
so large that the superconducting aspects of the model 
become irrelevant and one crosses over to the regime gov- 
erned predominantly by the physics of antiferromagnets. 

It should be emphasized that the phase diagram pre- 
sented here is based on the (inhomogeneous) MF treat- 
ment of a 2-dimensional model. As such, it is susceptible 
to the quantum and thermal fluctuations that tend to 
destroy long-range order. For instance, the ICAF phase 
in Figure | in a real material is more likely to manifest 
itself as incommensurate AF fluctuations, rather than a 
pure phase. However, the effects of the 3rd dimension 
and impurity pinning may stabilize the MF phases at a 
sufficiently low temperatures, revalidating the MF phase 
diagram. 

Within our model, we find that increasing on-site re- 
pulsion leads to a suppression of superconductivity. For 
larger U , the pure d-wave superconducting region of the 
phase diagram (SC) shrinks, and the width of the su- 
perconducting stripes in the SSC region decreases. For 
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example, changing U from At to 5i (keeping V fixed) 
leads to a five-fold reduction of the superconducting or- 
der parameter on stripes. This is in contrast with the 
homogeneous mean- field results 0, there the d-wave su- 
perconductivity is independent of the magnitude of U. 
On the other hand, large values of the on-site repulsion, 
U > At, in the Hubbard model lead to diagonal stripes 
p6[ . Applying our model to diagonal stripes (U = 5t, 
V = — 0.9t), we find that superconductivity vanishes 
completely at the optimal stripe filling. Perhaps it is 
not a coincidence, that the insulating cuprates do in fact 
show diagonal stripes, as opposed to the lattice-aligned 
stripes in the superconducting cuprates E^j. 

There is an analogy between the role played by doping 
in our model, and the strength of the electron-phonon 
coupling, A, in the McMillan criterion for the maximum 
achievable T c in conventional superconductors fl8| . In 
McMillan's picture, increasing A favors superconductiv- 
ity. However, increasing A too much induces a structural 
transition, and hence changes the reference ground state. 
In our model, increasing doping brings the superconduct- 
ing stripes closer together, and hence enhances the global 
T c . However, increasing doping too much causes a tran- 
sition to the uniform state, with a subsequent monoton- 
ically decreasing T c as a function of doping. 

The topology of the phase diagram we have proposed 
appears to be relevant for the superconducting cuprates, 
such as LSCO and YBCO. The same simple model (for 
other parameter values) can produce other topologies 
as well. A non-trivial topology, which may be real- 
ized in a material with a weaker attractive coupling, is 
when the superconducting region is fully contained in the 
AF/ICAF region. Depending on the parameters, such a 
material may be a striped superconductor in a certain 
range of dopings. Also, similar physics may occur in or- 
ganic charge-transfer salts that show AF, ICAF and 
superconductivity under pressure which controls inter- 
chain coupling. There, the intrinsically anisotropic cou- 
pling is also important for the stabilization of the meso- 
scopic inhomogeneities p0f . 

In conclusion, we find that a simple one-band model 
with on-site repulsion and nearest-neighbor attraction, 
in an appropriate range of parameters, can simultane- 
ously sustain both incommensurate antiferromagnetism 
and inhomogeneous superconductivity. Prompted by this 
finding and utilizing well-known antiferromagnetic and 
superconducting phase diagrams, we have constructed a 
generic phase diagram that captures many of the phases 



observed in the cuprate and organic superconductors. 
An experimental test of the Josephson-coupled super- 
conductor proposed here (see also [^i)) can be performed 
by measuring the in-plane Josephson plasmon resonance. 
Although our simple model appears to capture much of 
the observed rich physics, however it can be readily elab- 
orated (with multiple electron bands, 3-dimcnsionality, 
long-range interactions, lattice coupling, etc.) for more 
quantitative comparisons with specific materials. 
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